Magnetic Fluctuations of Filled Skutterudites 
Emerging in the Transition Region between Singlet and Triplet States 
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In order to clarify magnetic properties of filled skutterudites, we analyze the Anderson model 
including seven / orbitals hybridized with an a u conduction band using a numerical technique. For 
n=2 corresponding to Pr-based filled skutterudites, where n is the local /-electron number, even 
if the ground state is a singlet, there remain significant magnetic fluctuations from a triplet state 
with a small excitation energy. This result can be understood by the fact that /-electron states are 
clearly distinguished as itinerant and localized ones in the filled skutterudite structure. This picture 
also explains the complex results for /-electron magnetic susceptibility and entropy for n=l~13. 

PACS numbers: 74.25.Ha, 75.20.Hr, 74.70.Tx 
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The recent discovery of heavy fermion superconduc- 
tivity in PrOs4St>i2 Q has triggered intensive studies on 
filled skutterudite compounds in the research field of con- 
densed matter physics. Among the crystalline electric 
field (CEF) states of filled skutterudites characterized by 
the Th point group aTi singlet has been found exper- 
imentally to be the ground state of PrOs4Sbi2 |3j. How- 
ever, the results of NMR experiments Q suggest that the 
superconductivity is unconventional. In fact, PrOs4Sbi2 
exhibits exotic features such as multiple superconduct- 
ing phases 0] and the breaking of time-reversal symme- 
try as detected by /iSR experiments |(|. On the other 
hand, in the same family of Pr-based filled skutterudites, 
PrRu4Sbi2 is a conventional s-wave superconductor Q, 
PrOs4Pi2 is a non-magnetic metal ||, PrRu4Pi2 exhibits 
a metal-insulator transition ||, and PrFe4Pi2 shows ex- 
otic quadrupolar ordering . It is quite interesting that 
the electronic properties are easily changed by simple 
substitution of transition metal atoms and pnictogens. 

Besides Pr-based filled skutterudites, many other kinds 
of filled skutterudite compounds have been also synthe- 
sized. Those materials also exhibit a variety of elec- 
tronic properties: La-based filled skutterudite materials 
are known to be conventional BCS superconductors. Ce- 
bascd filled skutterudites are Kondo semiconductors with 
energy gaps up to a thousand Kelvins. For skutterudites 
containing rare-earth ions other than La, Ce, and Pr, a 
ferromagnetic ground state has been frequently observed, 
as for instance in RFe4Pi2 with R = Nd, Sm, Eu, Gd, Tb, 
Dy, and Ho. However, antiferromagnetic ground states 
have also been found in GdRu4Pi2 and TbRu4Pi2 [l0(. It 
is clearly desirable to identify the key issues which control 
the electronic properties of filled skutterudite compounds 
from a microscopic point of view. 

However, a microscopic theory for magnetism and su- 
perconductivity of /"-electron systems with n>l (n is 
local /-electron number) has not been satisfactorily de- 
veloped up to now, owing to the complexity of the many- 
body problem which stems from the competition among 
strong spin-orbit coupling, Coulomb interactions, and the 



CEF effect. To include such interactions, the LS coupling 
scheme has been widely used. However, that approach 
cannot exploit standard quantum-field theoretical tech- 
niques, since Wick's theorem does not hold in the LS 
coupling scheme. In order to overcome such a difficulty, 
we have proposed to construct a microscopic model for 
/-electron systems based on a j-j coupling scheme |l II . 
When we attempt to apply the latter model to filled skut- 
terudites, there are several issues to be addressed. In 
particular, one may have doubts about the application 
of the j-j coupling scheme to rare-earth materials, since 
Coulomb interactions are generally larger than the spin- 
orbit coupling in 4/-electron systems. It is also necessary 
to clarify how the CEF energy levels for /"-electron sys- 
tems are reproduced in the j-j coupling scheme. 

In this Letter, we analyze a multi-orbital Ander- 
son model, which correctly contains spin-orbit coupling, 
Coulomb interactions, and the CEF effect. For n=2, sig- 
nificant magnetic fluctuations are found to remain even 
in an / 2 -electron system with a Ti singlet ground state, 

(2) 

since there is a excited state triplet with a small ex- 
citation energy controlled by the CEF interaction. It is 
also shown that essentially the same result is obtained 
in an Anderson model constructed with the j-j coupling 
scheme. We find that in filled skutterudites, a certain / 
orbital hybridizes with the conduction band, while other 
orbitals are localized. This picture also explains the mag- 
netic properties of filled skutterudites for n=l~13. 
The Anderson Hamiltonian is written as 

H = V] £fe4 CT c kCT + V (V m clj mu + h.c.) + Hi, (1) 
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where is the conduction electron energy dispersion, 
Cka is the annihilation operator for conduction electrons 
with momentum k and spin a, a=+l (—1) for up (down) 
spin, f ma is the annihilation operator for /-electrons with 
spin a and angular momentum m(=— 3,- • -,3), and V m is 
the hybridization between / electrons and the conduc- 
tion band. In filled skutterudites, the conduction band is 
given by a u , constructed from p-orbitals of the pnictogens 
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FIG. 1: (a) Eigen-energies of Hf vs. x for n—2. (b) Excitation 
energy vs. x around at x~0.37. 

[l2| . Note that hybridization occurs between states with 
the same symmetry. Since the a u conduction band has 
xyz symmetry, we set V2=—V—2=V and zero for the other 
values of m. Throughout this paper, the energy unit is 
half of the bandwidth of the conduction band [l3( • 
The /-electron term Hf in Eq. is given by 

-^f ^ ^ (-^m,m' &aa' "h ^C?7l, £7,771' ,<j' ^) frnafm' 

+ ^m 1 ,m2,m 3 ,m i fm 1< j 1 fm 2rT2 frn 3 (T2fm 4 a 1 , (2) 
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where B m _ m i is the CEF potential for angular momentum 
I = 3, <S ffff ' = (1 + (T(t')/2, A is the spin-orbit coupling. 
The matrix elements are given by Cm,a,m,a = ma/2, 
Cm+i,Ti-m.±i = m(m ± l)/2, and zero for the 

other cases. The Coulomb integral /mi,m 2 , 777.3, m 4 is ex- 
pressed by the combination of four Racah parameters, 
A, £?, C, and D [14j . For the Th point group, -B m , m ' 
is given by three parameters, a;, y, and VF 0, where x 
and y specify the CEF scheme and W denotes the energy 
scale for the CEF potential. 

Let us first examine the local /-electron states, in or- 
der to set the parameters and check the validity of Hf. In 
Fig. 1(a), eigen-energies of Hf in the low-energy region 
are plotted as functions of x for n=2. The other CEF 
parameters, Racah parameters, and the spin-orbit inter- 
action are appropriately chosen for filled skuttcrudites. It 
is observed that the results for the LS coupling scheme 
are well reproduced, but we also emphasize that Hf cor- 
rectly reproduces the local /-electron states for both the 
LS and j-j coupling schemes, depending on the Coulomb 
interactions and spin-orbit coupling, for any value of n. 

(2) 

Recall that the Ti singlet ground state and T\ ' triplet 
excited state are split by only a small energy difference in 
Pr-based filled skutterudites ||- To study this, we mag- 

(2) 

nify the region in which the Tf and T 4 states cross. In 
the following, we consider the region in Fig. 1(b) where 
0.36 < x < 0.37. 

Here we include the hybridization between / electrons 
and conduction band states. For this purpose, we employ 
the numerical renormalization group (NRG) method [15| , 
in which the momentum space of the conduction electrons 
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FIG. 2: (a) Magnetic susceptibility and (b) entropy of / elec- 
tron vs. T for n=2 in the Anderson model Q for x=0.36, 
0.365, 0.366, and 0.37. (c) Critical values of x between mag- 
netic and non-magnetic regions for V=0.025, 0.05, 0.075 and 
0.1. Broken curve is drawn to fit the numerical results, (d) 
Eigen-energies of Hf for n=2. Insets denote /-electron con- 
figurations for Ei, T4, and Ts states, in which two / electrons 
are accommodated in the / 1 -energy levels, (e) Magnetic sus- 
ceptibility and (f) entropy of / electron vs. T for n=2 in the 
Anderson model @ for various values of A. 



is logarithmically discretized near the Fermi energy. The 
discretization is determined by a cut-off parameter A. In 
this paper we use A=5. Due to the limitation of com- 
puter memory, we keep only 4000 low-energy states for 
each renormalization step, but the results do not change 
qualitatively for A>3 when 4000 states are kept. 

In order to examine the properties of the / electrons, 
we evaluate the magnetic susceptibility \ and the entropy 
S of the / electrons. In Fig. 2(a), we show calculated re- 
sults for Tx for n=2 and several values of x between 0.36 
and 0.37. For x < 0.365, Tx goes to zero for small T, 
while for x > 0.366, Tx becomes constant at low temper- 
atures 0| . As shown in Fig. 2(b) for S, we obtain log 3 as 
residual entropy for £=0.366 and 0.37, indicating that a 
local triplet remains in the low-temperature region. Here 
we emphasize that magnetic fluctuations still remain at 
low temperatures, even when Ti is the ground state, if 
is an excited state with a small excitation energy. 
We suggest that anomalous properties of Pr-based filled 



3 



skutterudites originate from such magnetic fluctuations. 

We remark that the transition between magnetic and 
non-magnetic states is governed by the exchange inter- 
action J c f between / electrons and the conduction band, 
expressed as J c f=V 2 /SE, where SE denotes the energy 
difference between the f 2 and / 3 (or J 1 ) lowest-energy 
states. As shown in Fig. 2(c), the boundary curve be- 
tween magnetic and non-magnetic phases is proportional 
to V 2 , suggesting that the dominant energy scale should 
be J c f for the appearance of magnetic fluctuations. 

Note that the transition between magnetic and non- 
magnetic states seems to be abrupt, not gradual as ob- 
served in the usual Kondo system and the two-channel 
Anderson model [l7j . In the Kondo problem, the local 
moment is suppressed only by hybridization with conduc- 
tion electrons, but in the present case, a singlet ground 
state is also obtained through local level crossing due to 
the CEF potential, in addition to the hybridization pro- 
cess. Furthermore, localized orbitals exist in the present 
model, as will be discussed later. Namely, the magnetic 
moments of the / orbitals hybridized with the conduc- 
tion band are suppressed as in the Kondo effect, while the 
moments of the localized orbitals are not screened [lif . 
An abrupt change in the magnetic properties is caused 
by the duality of the / orbitals in combination with the 
level crossing effect due to the CEF potential. 

In order to understand the above point more explic- 
itly, it is useful to express the /-electron state with xyz 
symmetry in the j-j coupling scheme as 

|xyz, a) = V477I7/2, f 7 , or) - s/3/7\5/2, f 7 , a), (3) 

where \j,T,a) denotes the state in the j-j coupling 
scheme with total angular momentum 7, irreducible rep- 
resentation F for the Oh point group [19| , and pseudospin 
a [2p|. This relation indicates that only the f 7 state is 
hybridized with a u conduction band states with xyz sym- 
metry, leading to the Kondo effect, while the Tg electrons 
are localized and become the source of local fluctuations 
at low temperatures. It is an important issue of filled 
skutterudite structures that the nature of the / electrons 
can clearly be distinguished as itinerant r 7 or localized 
Tg. This point provides a possible explanation for the 
heavy-fermion phenomenon occurring in / 2 -electron sys- 
tems such as Pr-based filled skutterudites. 

In the j-j coupling scheme, the model is rewritten as 

H = E ^«L«k 5 + X> fl L/- + h - c -) + Hi, (4) 

where Ej- is the dispersion of the r 7 conduction electrons, 
ak& is the annihilation operator for conduction electrons 
with momentum k and pseudospin <r, / 7 a- is the anni- 
hilation operator for / electrons with pseudospin a and 
orbital 7 for the j=5/2 sextet on the impurity site, and V 
is the hybridization between conduction and / electrons 
with r 7 symmetry. The orbital index 7 distinguishes 



three kinds of Kramers doublets, two Tg (a and 6) and 
one r 7 (c) . The local /-electron term Hi is given by 

& = + E E C£ 2 ,73,7, 

7,<T 7l~74<7l,<T2 

X ./"ViCTl /72O2 f~t3&2 (5) 

where B 1 is the CEF potential and / is the Coulomb 
interaction in the j-j coupling scheme, expressed by other 
Racah parameters E , Ei , and E 2 [llj . 

Since the CEF potential is already diagonalized, it is 
convenient to introduce a level splitting A between P 7 
and Tg as A=£? r — Bp . In Fig. 2(d), we depict the low- 
energy states of H[ vs. A for n—2. The ground state for 
A<0 is a T 5 triplet composed of a couple of Tg electrons, 
while for A>0, it is a Ti singlet which is mainly composed 
of two P 7 electrons. Note that for A>0, the first excited 
state is a T4 triplet composed of r 7 and Tg electrons. 

Now we again employ the NRG method to investigate 
the magnetic properties of the Anderson model in the 
j-j coupling scheme. In Figs. 2(e) and (f), we depict 
T\ and S for several values of A. Racah parameters are 
set as appropriate values for rare-earth compounds. It is 
observed that T\ becomes constant and residual entropy 
log 3 remains at low temperatures for A < 5 x 10 -5 , while 
Tx goes to zero rapidly for A > 6 x 10~ 5 . Note that the 
transition is also dominated by J c f. 

Let us consider the reason why the essential magnetic 
features of Pr-based filled skutterudites can be captured 
as described in the j-j coupling scheme. As shown in 
Fig. 2(d), when the ground state is a l\ singlet, there 
are two triplet excited states, T4 and T§. The T$ triplet 
is composed of a pair of Tg electrons, while the P4 triplet 
is composed of one P 7 and one Tg electron. Since Pg 
electrons do not hybridize with T 7 conduction electrons, 
the P5 triplet can survive. Note that T4 triplets in Th are 
given by the mixtures of T4 and in Oh. Such a mixing 
is not included in the j-j coupling scheme, but there is a 
P5 excited state. Thus, the local triplet still remains, as 
long as the excitation energy is smaller than J c f. 

Thus far, we have focused on the case of n=2, but 
we can further study the cases of n=l~13 based on the 
original Anderson model Q. In the point-charge pic- 
ture, the CEF parameters do not depend on n, since the 
CEF effect is due to electrostatic potential from ligand 
ions surrounding the rare earth. In the following, we fix 
x=0.366, in which the local f 2 ground state is a Ti sin- 
glet, but for which we have found significant magnetic 
fluctuations at low temperatures. In Figs. 3(a) and (b), 
the magnetic susceptibility and entropy for each case of 
n<7 are shown. For n=2, 4, and 6, the ground state is a 
Ti singlet. Except for the case of n=2, the excitation en- 
ergy is large and both magnetic and orbital fluctuations 
should be rapidly suppressed with decreasing tempera- 
ture. Thus, Tx immediately becomes zero at low tem- 
peratures for n=4 and 6. On the other hand, for n=l, 
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FIG. 3: (a) Magnetic susceptibility and (b) entropy of / elec- 
trons vs. T for n=l~6. (c) Magnetic susceptibility and (d) 
entropy for n— 7~13. We set 1=0.366 and all other parame- 
ters are the same as those in Figs. 2 (a) and (b). 



r§ states are localized. For n=2, corresponding to Pr- 
based filled skutterudites, it has been found that mag- 

(2) 

netic fluctuations significantly persist, if there is a I4 
excited state triplet with a small excitation energy. It 
has been shown that essentially these same results are 
reproduced by the Anderson model in the j-j coupling 
scheme. We believe that these results open a new door 
for the study of magnetic and superconducting proper- 
ties of Pr-based filled skutterudites from a microscopic 
point of view. 
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cility of the Supercomputer Center of Institute for Solid 
State Physics, University of Tokyo. 



3, and 5, the ground state is a Tq7 quartet in Th (=r§ in 
Oh), as confirmed from the residual entropy of log 4. 

In Figs. 3(c) and (d), we show the NRG results for sus- 
ceptibility and entropy for n>7. First of all, the absolute 
values of x are much larger than those for n<7, since the 
total angular momentum J becomes large for n>7 owing 
to Hund's rule coupling. Typically, at half-filling, total 
spin S(—J) is equal to 7/2, and the Curie constant for an 
isolated ion is as large as 21/ig/fcB- Over a broad tem- 
perature region, this value has been observed in Fig. 3(c) 
for n=7, indicating that the S=7/2 spin survives at rel- 
atively low temperatures. In fact, we clearly see an en- 
tropy of log 8 from the 5=7/2 octet, as shown in Fig. 3(d) 
[21|. For the cases of n=8 and 12, the ground state is a 
Ti singlet and the magnetic excited state energy is now 
large, in sharp contrast with the case of n=2. Thus, the 
susceptibility rapidly goes to zero. For n=9, 11, and 13, 
the local ground state is a T 5 doublet in Th (=F 6 in Oh). 
Thus, this state does not hybridize with the conduction 
band and the magnetic moment still persists even in the 
low-temperature region. In fact, we observe a residual 
entropy of log 2 in these cases, as shown in Fig. 3(d). For 
n=10, the local ground state is a r 4 triplet, but as ob- 
served in Figs. 3(c) and (d), the local triplet seems to 
remain at low temperatures. This is easily understood, 
if we recall that T 4 triplets of Th & r e given by mixtures 
of T4 and r 5 of Oh- As mentioned above, the T 5 triplet 
still persists even after hybridization, since a r 8 electron 
does not hybridize with a conduction electron. 

In summary, we have analyzed a multi-orbital Ander- 
son model using the NRG method. We have established 
that Ft states hybridize with the conduction band, while 
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